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We present some arguments suggesting that the mismatch between Bekenstein- Hawking entropy 
and the entropy of entanglement for vector fields is due to the same gauge configurations which 
saturate the contact term in topological susceptibility in QCD. In both cases the extra term with a 
"wrong sign" is due to distinct topological sectors in gauge theories. This extra term has non-dispersive 
nature, can not be restored from conventional spectral function through dispersion relations, and can 
not be associated with any physical propagating degrees of freedom. We make few comments on some 
profound consequences of our findings. In particular, we speculate that the source of the observed 
dark energy may also be related to the same type of gauge configurations which are responsible for 
the mismatch between black hole entropy and the entropy of entanglement in the presence of causal 
horizon. 



I. INTRODUCTION 

The relation between black hole entropy and entropy 
of entanglement for matter fields has been a subject of in- 
tense discussions for the last couple of years, see reviews[l- 
4] and original references therein. There are many sub- 
tleties in relating these two things. The present work is 
concentrated just on one specific subtlety first discussed 
in [5]. Namely, it has been claimed [5] that for spins 
zero and one-half fields the one loop correction to the 
black hole entropy is equal to the entropy of entanglement 
while for spin one field the black hole entropy has an extra 
term describing the contact interaction with the horizon. 
Precisely this contact interaction with the horizon is the 
main topic of the present work. Before we elaborate on 
this subject we want to make one preliminary remark 
regarding the term "black hole entropy" . As it has been 
argued in a number of papers, see e.g. [6-8] and also review 
papers [1-4] the notion of black hole entropy should apply 
not just to black holes but to any causal horizon ("black 
hole entropy without black holes"). We adopt this view- 
point, and in fact we shall not discuss black hole physics 
in this paper at all. Rather, the main application of our 
studies will be cosmology of the expanding universe and 
its causal horizon. For short, we shall use term "entropy" 
through the paper. 

The unique features of the contact term related to the 
vector gauge field in the entropy computations can be 
summarized as follows [5] (see also follow up paper [9]) : 

a) the contact term being a total derivative can be 
represented as a surface term determined by the 
behaviour at the horizon; 

b) this term makes a negative contribution to the black 
hole entropy. 

c) therefore, it can not be identified with entropy of 
entanglement which is intrinsically positive quantity. 

d) this contribution does not vanish even in two dimen- 
sions when the entropy of entanglement is identically 
zero as no physical propagating degrees of freedom 
are present in the system; 



e) this contribution is gauge invariant in two dimen- 
sions and gauge dependent in four dimensional case 

[9]; 

f ) the technical reason for this phenomenon to happen 
is as follows. One can not use the physical Coulomb 
gauge (when only physical degrees of freedom are 
present in the system) as it breaks down at the 
origin, where Ag is ill defined. Therefore, an al- 
ternative description in terms of a covariant gauge 
(when unphysical degrees of freedom inevitably ap- 
pear in the system) should be used instead. 

g) in this covariant description the entropy is obtained 
by varying the path integral with respect to the 
deficit angle of the cone as explained in [5] . Such 
a procedure can (in principle) lead to a negative 
value for the entropy. In fact, it does come out 
negative [5]. 

The main goal of the present work is to argue that the 
presence of this "weird" term is intimately related to the 
well known property of gauge theories when a summation 
over all topological sectors must be performed for the path 
integral to be properly defined. We explain how all fea- 
tures a)-g) from the list above can be natural understood 
within our framework when summation over topological 
sectors is properly taken into account. In Minkowski 
space the corresponding procedure is known to produce a 
non-dispersive contact term with a "wrong sign" which 
plays a crucial role in resolution of the so-called U(1)a 
problem in QCD. Exactly this non-dispersive contact term 
eventually becomes the "weird" term with properties a)- 
g) when we go from conventional Minkowski space into 
curved/time-dependent background with a causal horizon. 

Our consideration in this paper will be based on analysis 
of the local characteristics (such as topological suscepti- 
bility, free energy density, etc) computed deep inside the 
horizon region. It is very different from the computation 
of a global characteristic such as the total entropy for a 
black hole when the closest vicinity of the horizon (just 
outside of it) plays the crucial role in the computations. 
Nevertheless, we remain sensitive to the existence of hori- 
zon because our analysis is based on consideration of some 
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specific topologically- protected quantities. Essentially, 
by analyzing very unusual features a)-g) listed above we 
learn some important lessons on behaviour of the ground 
state resulting from merely existence of a causal horizon 
in the presence of the gauge degrees of freedom in the 
system. 

The workflow is as follows. In section II we present 
our arguments for two dimensional case when all compu- 
tations can be explicitly performed. We generalize our 
arguments for four dimensional case in section III. We 
argue that this term is indeed gauge dependent in four 
dimensions in abelian case as explicit computations of 
ref. [9] suggest. However, we shall argue that this term 
becomes a gauge independent in non-abelian case. We 
make few comments on some profound consequences of 
our findings in section IV, where we speculate that the 
source of the observed dark energy might be related to 
the same gauge configurations which are responsible for 
the mismatch between black hole entropy and the entropy 
of entanglement. 



II. TOPOLOGICAL SECTORS, CONTACT 
TERM WITH "WRONG SIGN", AND ALL THAT 
FOR 2D QED IN RINDLER SPACE 

First of all, we shall demonstrate below the presence 
of a nonconventional contribution into the energy with 
a "wrong sign" in Minkowski space. This contribution 
is gauge invariant, it exists even in pure photo-dynamics 
when no propagating degrees of freedom are present in 
the system. It can not be removed by any means (such 
as redefinition of the energy) as it is a real physical con- 
tribution. In particular, the anomalous Ward Identities 
(which emerge when the massless fermions are added into 
the system) can not be satisfied without this term. We 
shall argue that this term can be treated as a contact 
term, and in fact is related to the existence of different 
topological sectors in this (naively trivial) two dimen- 
sional photo-dynamics. In different words, the presence of 
different topological sectors in the system, which we call 
the "degeneracy" for short 1 , is the source for this contact 
term which is not related to any physical propagating 
degrees of freedom. 

As the next step, we shall discuss the same system in 
the presence of the horizon in the Rindler space. We shall 



Not to be confused with conventional term "degeneracy" when two 
or more physically distinct states are present in the system. In the 
context of this paper the "degeneracy" implies there existence of 
winding states \n) constructed as follows: T\n) = \n + 1). In this 
formula the operator T is the large gauge transformation operator 
which commutes with the Hamiltonian [T, H] = 0. The physical 
vacuum state is unique and constructed as a superposition of 
|n) states. In path integral approach the presence of n different 
sectors in the system is reflected by summation over k £ Z in eq. 
(10, 11,12). 



argue that the contact term (which emerges itself as a re- 
sult of topological features of gauge theory) demonstrates 
the "weird" and strange properties listed above in the 
presence of the horizon. 

In what follows it is convenient to study the topological 
susceptibility x (rather than free energy itself) which is 
related to the 9 dependent portion of the free energy 
density 2 as follows 

fa a n\ <9 2 F vac (/3, 0) 
x(/3,6» = 0) = — , (1) 

uu 9=0 

where 9 is the conventional 9 parameter which enters 
the Lagrangian along with topological density operator, 
see precise definition below. We always assume that 
9 = 0, however x(& — 0) ^ does not vanish, and in 
fact is the main ingredient of the resolution of the U(1)a 
problem in QCD [10-12], see also [13-15]. Free energy 
itself F vac (9) can be always restored from \ as dependence 
on 9 is known to be -Fvac ~ cos 9. As we show below, the 
topological susceptibility \ (and therefore F vac ), being 
the local characteristics of the system, nevertheless are 
quite sensitive to merely existence of horizon, even when 
computed far away from it. As we shall see this sensitivity 
is related to the degeneracy of the system and topological 
nature of \- 

A. Topological susceptibility and contact term 

The simplest (and physically attractive) choice is 
Coulomb gauge when no physical propagating degrees 
of freedom are present in the system, and therefore the 
dynamics must be trivial. It is well known, why this naive 
argument fails: the vacuum in this system is degener- 
ate, and one should consider an infinite superposition of 
of the winding states |n) as originally discussed in [16]. 
Such a construction in Coulomb gauge restores the cluster 
and other important properties of quantum field theory. 
The vacuum in this gauge is characterized by long range 
forces (if charged physical fermions are introduced into 
the system). This long range force prevents distant re- 
gions from acting independently. We believe that precisely 
this feature leads to the difficulties mentioned in [5] in 
computations of the entropy in physical Coulomb gauge 
in two dimensions, where a covariant gauge has been used 
instead. 

As our goal is to make a connection with computations 
of ref. [5] we shall not elaborate on Coulomb gauge in the 
present paper any further, but rather consider a covariant 
gauge to study this system. In the covariant Lorentz 
gauge there are no long range forces. Instead, a new 



2 in case of infinite manifold (rather than finite size /3 = T _1 ) the 
free energy from relation (1) becomes the conventional vacuum 
energy as employed in study of the U(1)a problem in QCD in 
[10-12]. 
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(unphysical) degrees of freedom emerge in the system, see 
precise definition below. 

We want to study the topological susceptibility x m 
the Lorentz gauge defined as follows 3 , 



X 



Air 2 fc-s-o 



lim / d 2 xe lkx (TE{x)E{0)) , 



(2) 



where Q — ^E is the topological charge density and 



d 2 x Q{x) 



e 

2^ 



d 2 x E(x) = k 



(3) 



is the integer valued topological charge in the 2d U(l) 
gauge theory, E(x) — d\A 2 — d%A\ is the field strength. 
The expression for the topological susceptibility in 2d 
Schwinger QED model is known exactly [17] 



Xqed 



4tt 2 



d 2 x 



S 2 (x) 



2tt 2 



K (n\x\) 



, (4) 



where /x 2 = e 2 /ir is the mass of the single physical state in 
this model, and Kq([i\x\) is the modified Bessel function 
of order 0, which is the Green's function of this massive 
particle. The expression for \ for pure photo-dynamics is 
given by (4) with coupling e = in the brackets 4 which 
corresponds to the de-coupling from matter field ip, i.e. 



Xeum 



4tt 2 



d 2 a; 



[S 2 ( a 



(5) 



The crucial observation here is as follows: any physical 
state contributes to \ with negative sign 



Xdispersive 



HmV 



-E\n)(n\^E\0) 
—k 2 — m 2 



<0. (6) 



In particular, the term proportional —Ko(p,\x\) with neg- 
ative sign in eq. (4) is resulted from the only physical 
field of mass /j,. However, there is also a contact term 
/ d 2 x [<5 2 (x)] in eqs. (4), (5) which contributes to the 
topological susceptibility \ with the opposite sign, and 
which can not be identified according to (6) with any 
contribution from any physical asymptotic state. 

This term has fundamentally different, non-dispersive 
nature. In fact it is ultimately related to different topo- 
logical sectors of the theory and the degeneracy of the 
ground state [18] as we shortly review below. Without 
this contribution it would be impossible to satisfy the 
Ward Identity (WI) because the physical propagating de- 
grees of freedom can only contribute with sign (— ) to the 
correlation function as eq. (6) suggests, while WI requires 
X = in the chiral limit m = 0. One can explicitly check 



that WI is indeed automatically satisfied 5 only as a re- 
sult of exact cancellation between conventional dispersive 
term with sign (— ) and non-dispersive term (5) with sign 

(+), 



e 

4^ 
4tt 2 



d 2 x 



5 2 {x) 



■^K Q {n\x\) 



(7) 



1 - 



7T /i z 



47T 2 



[1-1]=0. 



B. 



The origin of the contact term— summation over 
topological sectors. 



The goal here is to demonstrate that the contact term 
in exact formulae (4), (5) is a result of the summation over 
different topological sectors in the 2d pure U{1) gauge 
theory as we now show. We follow [17] and introduce 
the classical "instanton potential" in order to describe 
the different topological sectors of the theory which are 
classified by integer number k defined in eq. (3). The 
corresponding configurations in the Lorentz gauge on two 
dimensional Euclidean torus with total area V can be 
described as follows [171: 



A (k) 



irk 



OEM 



2irk 



eV ^ 

such that the action of this classical configuration is 



d 2 xE 2 



2 1,2 



2lT 2 k 

e 2 V 



(8) 



(9) 



This configuration corresponds to the topological charge 
k as defined by (3). The next step is to compute the 
topological susceptibility for the theory defined by the 
following partition function 



Z = E / VAe-if* 2 * 1 



(10) 



All integrals in this partition function are gaussian and 
can be easily evaluated using the technique developed in 
[17]. The result is determined essentially by the classical 
configurations (8), (9) as real propagating degrees of 
freedom are not present in the system of pure t/(l) gauge 
field theory in two dimensions. We are interested in 
computing \ defined by eq. (2). In path integral approach 
it can be represented as follows, 



4tt 2 Z 



d 2 xE(x)E(0)e 



-\JSxE 2 



(11) 



3 Here we use Euclidean metric where path integral computations 
(4) have been performed. 

A 2 

factor in front of (4) does not vanish in this limit as it is due 
to our definition (2) rather than result of dynamics 



5 When m ^ the WI takes the form \ ~ m\{ipip)\. It is also 
automatically satisfied because fi 2 = — +0(m), and cancellation 
in eq. (7) is not exact resulting in behaviour \ ~ m m complete 
accord with WI [19]. 
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This gaussian integral can be easily evaluated 6 and the 
result can be represented as follows [18], 



4tt 2 



• V ■ 



4n 2 k 2 
e 2 V 2 



exp(- 



2n 2 k 2 \ 
e 2 V I 



Efeez ex P(" 



2ir 2 k 2 
e 2 V 



(12) 



In the large volume limit V — > oo one can evaluate the 
sums entering (12) by replacing J2 keZ ~^ J dk such that 
the leading term in eq. (12) takes the form 



X 



e 

4^2 



V ■ 



4tt 2 
e 2 V 2 



e 2 V 
4^2 



e 

4^2' 



(13) 



Few comments are in order. First, the obtained expression 
for the topological susceptibility (13) is finite in the limit 
V — y oo , coincides with the contact term from exact 
computations (4), (5) performed for 2d Schwinger model 
in ref. [17] and has "wrong" sign in comparison with 
any physical contributions (6). Second, the topological 
sectors with very large k ~ \J e 2 V saturate the series 
(12). As one can see from the computations presented 
above, the final result (13) is sensitive to the boundaries, 
infrared regularization, and many other aspects which are 
normally ignored when a theory from the very beginning is 
formulated in infinite space with conventional assumption 
about trivial behaviour at infinity. Last, but not least: the 
contribution (13) does not vanish in a trivial model when 
no any propagating degrees of freedom are present in the 
system! This term is entirely determined by the behaviour 
at the boundary, which is conveniently represented by the 
classical topological configurations (8) describing different 
topological sectors (3), and accounts for the degeneracy 
of the ground state 7 . We know that this term must be 
present in the theory when the dynamical quarks are 
introduced into the system. Indeed, it plays a crucial role 
in this case as it saturates the WI as formula (7) shows. 



C. The ghost as a tool to describe the contact term 

The goal here is to describe exactly the same contact 
term (5, 13) without explicit summation over different 
topological sectors, but rather, using the auxiliary ghost 
fields as it has been originally done in ref. [16] (using 
the so-called the Kogut-Susskind dipole). This auxiliary 
ghost field effectively accounts for the degeneracy of the 
ground state as discussed above. The computations in 
both refs. [5] and [16] are performed precisely in terms 
of the same auxiliary scalar field defined as follows 



(14) 



6 One can check that the contribution resulting from the quantum 
fluctuations about the classical background (8) does not change 
the result (13). Indeed, the corresponding extra "quantum" con- 
tribution ~% ■ J d 2 x[8 2 (x) — y] = vanishes as expected. 

7 See footnote 1 for clarification of the term "degeneracy" . In the 
given context the degeneracy implies the summation over k £ Z 
in eq. (10, 11,12). 



This formal connection allows us to make a link between 
expressions (5, 13) for the contact term with "wrong" sign 
computed in our framework in terms of the auxiliary scalar 
field as described below and the entropy computations 
performed in ref. [5] featuring the "weird properties" e)-g) 
as listed in Introduction. 

Our starting point is the effective Lagrangian describ- 
ing the same two dimensional gauge system. However, 
now the theory is formulated in covariant Lorentz gauge 
in terms of the scalar fields [16]. The crucial element 
accounting for different topological sectors of the underly- 
ing theory, and corresponding degeneracy of the ground 
state, does not go away in this description. Rather, this 
information is now coded in terms of unphysical ghost 
scalar field which provides the required "wrong" sign for 
contact term (5, 13). 

Precise construction goes as follows. The effective La- 
grangian describing the low energy physics (in Minkowski 
metric) is given by [16]: 



(15) 



1 



\i 4> + m\(qq) | cos 2y/TT 



The fields appearing in this Lagrangian are 



the ghost, 



its partner , 



(16) 



while <f> is the only physical massive degree of freedom. It 
is important to realize that the ghost field <j)\ is always 
paired up with <f>2 in each and every gauge invariant matrix 
element, as explained in [16]. The condition that enforces 
this statement is the Gupta-Bleuler-like condition on the 
physical Hilbert space H p hy S which reads like 



(h - 4>i) i+) \n phys ) = o, 



(17) 



where the (+) stands for the positive frequency Fourier 
components of the quantized fields. One can easily un- 
derstand the origin for a wrong sign for the kinetic term 
for (fii field. It occurs as a result of D 2 operator when the 
Maxwell term E 2 ~ D 2 is expressed in terms of the scalar 
field (14). As usual, the presence of 4-th order operator is 
a signal that the ghost is present in the system. Indeed, 
the relevant operator [□□ + /z 2 D] which emerges for this 
system can be represented as the combination of the ghost 
4>\ and a massive physical <j> using the standard trick by 
writing the inverse operator as follows 



□□ + n 2 a 



-□ 



1 



(18) 



This is a simplified explanation on how sign (— ) emerges in 
Lagrangian (15) describing auxiliary <pi field, see original 
paper[16] for details. 

The contact term in this framework is precisely rep- 
resented by the ghost contribution[18, 19] replacing the 
standard procedure of summation over different topologi- 
cal sectors as discussed above II B. Indeed, the topolog- 
ical density Q = ^-E in 2d QED is given by ^-E — 



5 



(^)^T ~ D ^i) [ 16 1- The relevant correlation func- 
tion in coordinate space which enters the expression for 
the topological susceptibility (2) can be explicitly com- 
puted using the ghost as follows 



Z7T Z7T 



f-V 

1 



TV 



d 2 P 



(27T) 



2 p 4 e- lpx 



5 2 (x) 



e 

2t^ 



K (p\x\) 



(19) 



where we used the known expressions for the Green's 
functions (the physical massive field ij> as well as the ghost 
4>i field) determined by Lagrangian (15) and switched 
back to Euclidean metric for comparison with previous 
results from sections II A, II B. 

The obtained expression precisely reproduces exact 
result (4) as claimed. In the limit e — > when matter fields 
decouple from gauge degrees of freedom we reproduce the 
contact term (5, 13) which was previously derived as a 
result of summation over different topological sectors of 
the theory. The non-dispersive contribution manifests 
itself in this description in terms of unphysical ghost 
scalar field which provides the required "wrong" sign for 
contact term. 

At the same time, this unphysical ghost scalar field does 
not violate unitarity or any other important properties of 
the theory as consequence of Gupta-Bleuler-like condition 
on the physical Hilbert space (17). Indeed, while the 
ghost's number density operator, N may look dangerous 
due to the sign (— ) in the commutation relations 



N = 



h,bl\=6. 



J2 { b l bk - 4«fe) 



(20) 



= -8 



kk' 



one can in fact check that the expectation value for any 
physical state vanishes as a result of the subsidiary condi- 
tion [16, 18]: 

(H phys |N|H phys > = , (o fc - b k ) \H phys ) = 0. (21) 

This vanishing result (21) obviously implies that no en- 
tropy may be produced in Minkowski space. In different 
words, the fluctuations of unphysical fields described by 
operator (20) do not lead to any physical consequences 
(except for merely existence of the contact term (5) as 
already discussed). 

We shall see in next subsection how this simple picture 
drastically changes when we consider the very same sys- 
tem but in presence of the horizon. We shall argue that 
the number density N of "fictitious particles" with wrong 
commutation relations starts to fluctuate in the presence 
of the horizon, in contrast with eq. (21). Therefore, we 
formulate a conjecture that precisely these fluctuations 
are responsible for term with a "wrong sign" in entropy 
computations [5, 9]. The corresponding contribution, as 



we already mentioned, is not related to any physical prop- 
agating degrees of freedom but rather, is due to presence 
of topological sectors in gauge theories (and the degener- 
acy of the ground state as its consequence, see footnote 
1 for clarification on terminology) which eventually lead 
to a non-dispersive contribution in topological suscepti- 
bility. To simplify things in what follows we consider a 
simple Rindler space when the Bogolubov's coefficients 
are exactly known. However, we argue that a generic 
case (when horizon is present in the system) leads to very 
similar conclusion. 



D. Rindler space 

The total entropy with "weird" properties listed in 
Introduction was computed a while ago [5], and there 
is no reason to review these results in the present paper. 
These original results have been reproduced in [9] by 
using another technique. Furthermore, in the same paper 
[9] it has been demonstrated that in two dimensional 
case the final result is gauge invariant, and therefore, it 
obviously represents a physically observable quantity. As 
we mentioned earlier, we are not interested in computing 
the global characteristics such as total entropy. Rather, 
we are interested in computing some local properties, such 
as topological susceptibility, or 9— dependent portion of 
the energy density (1) in the presence of the horizon. 
However, we shall argue below, the source for "weird" 
features in both cases is the same, and, in fact, related to 
fundamental properties of gauge theories as discussed in 
section II B. 

As we explained above, the presence of different topo- 
logical sectors in gauge theory (and the degeneracy of the 
ground state as its consequence) leads to contact term (5) 
even when no physical propagating degrees of freedom are 
present in the system. In physical Coulomb gauge this 
term manifests itself as the presence of long range force 
which prevents distant regions from acting independently. 
The same feature but in covariant Lorentz gauge is ex- 
pressed in terms of new (unphysical) degrees of freedom 
(16) which emerge in the system and effectively repro- 
duce the contact term as explicit computations show (19). 
While these unphysical degrees of freedom do fluctuate, 
these fluctuations do not lead to any physical observable 
expectation values in Minkowski space (21) as a result 
of cancellation between two unphysical fields similar to 
conventional Gupta-Bleuler condition in QED when two 
unphysical photon's polarizations cancel each other. We 
want to see how this conclusion changes when a horizon 
is present in the system. 

One can repeat the construction of previous section 
II C to describe (unphysical) degrees of freedom but in 
Rindler space [18]. A Rindler observer in (R,L) wedge will 
measure the number density of unphysical states using 
density operator N( fl ' L ) which is given by 

NW)=^(^iW>- a W]t a W)) . (22) 



G 



The subsidiary condition (17) defines the physical sub- 
space for accelerating Rindler observer 



'(R,L) AR,L) 
. a k u k 



H 



(R,L) 
phys 



= o, 



(23) 



such that the exact cancellation between </>i and (f>2 fields 
holds for any physical state defined by eq. (23), i.e. 



H 



(R,L) 
phys 







(24) 



as it should. However, if the system is prepared as the 
Minkowski vacuum state |0) then a Rindler observer us- 
ing the same operator for N( R ' L ) (22) will observe the 
following number density in mode k, 



(0|N ( ^ L) |0) = (0| ( 



h {R,L)\ h {R,L) 



a[ R ^ai R ^) |0) 



2 e ~™/a 



(25) 



where we used known Bogolubov's coefficients mixing 
the positive and negative frequency modes for operators 
b^f" , ai 1 '^ describing unphysical fluctuations [18]. 

One can explicitly see why the cancellation (21) of un- 
physical degrees of freedom in Minkowski space fail to hold 
for the accelerating Rindler observer (25). The technical 
reason for this effect to occur is the property of Bogol- 
ubov's coefficients which mix the positive and negative 
frequencies modes. The corresponding mixture can not 
be avoided because the projections to positive -frequency 
modes with respect to Minkowski time t and positive 
-frequency modes with respect to the Rindler observer's 
proper time rj are not equivalent. The exact cancellation 
of unphysical degrees of freedom which is maintained 
in Minkowski space can not hold in the Rindler space 
because it would be not possible to separate positive fre- 
quency modes from negative frequency ones in the entire 
spacetime, in contrast with what happens in Minkowski 
space where the vector d/dt is a constant Killing vector, 
orthogonal to the t — const hypersurface. The Minkowski 
separation is maintained throughout the whole space as 
a consequence of Poincare invariance. It is in a drastic 
contrast with the accelerating Rindler space [18]. 

The nature of the effect is the same as the conven- 
tional Unruh effect [20] when the Minkowski vacuum |0) 
is restricted to the Rindler wedge with no access to the 
entire space time. An appropriate description in this case, 
as is known, should be formulated (for R observer) in 
terms of the density matrix by "tracing out" over the 
degrees of freedom associated with L-region. In this case 
the Minkowski vacuum |0) is obviously not a pure state 
but a mixed state with a horizon separating two wages, 
which is the source of the entropy. In contrast with Unruh 
effect [20], however, one can not speak about real radiation 
of real particles as the ghost <fii and its partner (j>2 are 
not the asymptotic states and the corresponding positive 
frequency Wightman Green function describing the dy- 
namics of these fields vanishes [18]. In different words, 



these auxiliary fields contribute to the non-dispersive por- 
tion of the correlation function in eqs. (4,5,7), but not 
to conventional dispersive part which is unambiguously 
determined by the absorptive function as conventional 
dispersion relation dictates. 

Few more comments on (25) are in order. The effect 
is obviously sensitive to the presence of the horizon, and, 
therefore is infrared (IR) in nature. The IR nature of the 
effect was anticipated from the very beginning as formu- 
lation of the problem in terms of auxiliary fields (16) is 
simply a convenient way to deal with different topologi- 
cal sectors of the gauge theory in covariant gauge (and 
the degeneracy of the ground state as their consequence) 
instead of dealing with the long range forces in the uni- 
tary Coulomb gauge as discussed in sections II B and 
IIC. Also: the contribution of higher frequency modes 
are exponentially suppressed ~ exp(— to/a) as expected. 
The interpretation of eq. (25) in terms of particles is very 
problematic (as usual for such kind of problems) as typical 
frequencies when the effect (25) is not exponentially small, 
are of order w ~ a, and notion of "particle" for such uj is 
not well defined. 

We do not attempt to reproduce known results on 
entropy from ref. [5] based on non-vanishing expectation 
value for number density operator (25). First of all, it is 
not obvious what would be the physical meaning of such 
a computation based on expectation value (25) for the 
operator which satisfies "wrong" commutation relation 
(20). Furthermore, it is not obvious how to interpret N 
particles from eq.(25) when entire notion of particles is 
not even defined for relevant parameters. Indeed, as we 
argued above the effect is large (0|N|0) ~ 1 only for very 
large wave length A > a -1 which is the size of the horizon 
scale. 

Our goal here is in fact quite different. We want to 
argue that the source for the "wrong sign" in entropy 
computations [5] (featuring the "weird properties" as 
listed in Introduction) and the source for the "wrong sign" 
for the contact non-dispersive term (discussed in present 
paper) are in fact have the same origin. In addition to the 
arguments presented above, we note that the technical 
computations of the entropy performed in [5] are actually 
based precisely on the same representation for field (14) 
describing fluctuations of unphysical auxiliary degrees of 
freedom. This representation for field in our formalism 
eventually leads to the expression for the contact non- 
dispersive contribution ~ 8 2 {x) with "wrong" sign (19) 
and non vanishing number density (25) while in ref. [5] the 
very same representation for A^ field (14) leads to the 
"wrong" sign for the entropy. Furthermore, the contact 
term (5) can be represented as a surface term, 



Xe&m 



,(26) 



analogous to the "weird" contribution in the entropy com- 
putations [5, 9]. It is important to realize that the contact 
term (5, 13, 26) is a result of summation over all topo- 
logical sectors with inclusion of all quantum fluctuations 
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which account for the degeneracy of the ground state as 
discussed in section II B. At the same time, quite mirac- 
ulously, the final result (5, 13, 26) can be interpreted as 
a surface integral of a single classical configuration of a 
pure gauge field A c ^ ~ d^cf) 01 defined on a distance surface 
S\ and characterized by unit winding number 

r A c Ml^ r rAfd.6 [ d9 d<j> cl 

i -v - i -t- -J s-s- (27) 

These observations strongly suggest that the term with 
a "wrong" sign in the expression for entropy derived in 
[5] has exactly the same origin as the "wrong" sign for 
the contact term (5) as in both cases the relevant physics 
is determined by the surface integrals, not related to any 
physical propagating dynamical degrees of freedom. Fur- 
thermore, in both cases the sign of the effect is opposite 
to what one should expect from physical degrees of free- 
dom, and, finally, in both cases the starting point (formal 
representation for A M field (14)) is the same. 

• Therefore, we conjecture that: the surface term with 
a "wrong sign" in entropy computations [5, 9] and the 
"wrong sign" in topological susceptibility (2,5,13) are both 
originated form the same physics, and both related to 
the same (topologically nontrivial) gauge configurations, 
and must be present (or absent) in both computations 
simultaneously. In both cases a "wrong sign" emerges due 
to unphysical degrees of freedom fluctuating in far infrared 
(IR) region. The technical treatments of these terms in our 
framework and in ref. [5] of course is very different: we use 
conventional Hamiltonian approach supplemented by the 
condition (17) while in ref. [5] the Rindler Hamiltonian is 
ill defined on the cone, and computations are performed 
using some alternative methods. Nevertheless, in our 
framework, we interpret the fluctuations (25) of "fictitious 
particles" with "wrong" commutation relations (20) as a 
different manifestation of the same physics which led to a 
wrong sign in entropy computations [5, 9]. An additional 
argument supporting our conjecture will be presented in 
the next section where we show that these very different 
quantities nevertheless behave very similarly when the 
system is generalized from two to four dimensions, and 
therefore, they must be originated from the same physics. 

Our final comment here is this: the IR physics pene- 
trates into the physical gauge invariant correlation func- 
tion (2) not due to the massless degrees of freedom in the 
physical spectrum (there are none in fact), but rather, as 
a result of degeneracy of the ground state and summation 
over all topological sectors in gauge theory as discussed 
in sections II A, II B. The ghost (16) in this framework 
is simply a convenient tool to account for this far IR 
physics as it effectively accounts for the non-dispersive 
contact term with a "wrong sign" (19). It fluctuates in 
the presence of the horizon (25), and responsible for a 
"wrong" sign in entropy computations, according to our 
conjecture. However, it remains unphysical auxiliary field 
as it does not belong to the physical Hilbert space (and it 
never becomes an asymptotic state capable to propagate 
to infinity) [18]. It is interesting to notice that there are 



other known examples when the degeneracy of the ground 
state in the presence of the horizon leads to mismatch 
between black hole entropy and entropy of entanglement, 
see Appendix A for references and details. 



III. GENERALIZATION TO 4-D CASE 

The goal of this section is twofold. First, in next sub- 
section III A we make few comments on generalization of 
two dimensional results discussed above to four dimen- 
sional QED. In this case the corresponding calculations 
of the entropy are known [5, 9]. Analysis of these results 
further support our conjecture on common nature of the 
surface term with a "wrong sign" in entropy computations 
[5, 9] and "wrong sign" in topological susceptibility as 
the behaviour of the system follows precisely the pattern 
dictated by the conjecture. Secondly, in section III B we 
discuss four dimensional non-abelian gauge theories when 
corresponding computations of the entropy are not yet 
known. Nevertheless, based on our conjecture on com- 
mon origin of these two different phenomena, we predict 
a possible outcome if the corresponding computations are 
performed. 



A. Four dimensional abelian QED 

We start by reviewing the basic results of refs. [5, 9] 
on entropy computations in four dimensional case. In 
original paper [5] the gauge invariance of the "surface 
term" has not been tested. This question has been specif- 
ically discussed in followup paper [9] where it has been 
demonstrated that in two dimensions the result is indeed 
gauge invariant and coincides with the original expression 
found in ref. [5]. However, a similar analysis in four di- 
mensions turned out to be much more subtle, see details 
in [9]. In particular, it has been found that this term is 
gauge dependent in four dimensional abelian case, and 
therefore, it was discarded [9]. 

How one can understand such puzzling behaviour of 
the system when one jumps from two to four dimen- 
sions? If one accepts our conjecture formulated above 
then this puzzle has a very natural explanation. Indeed, 
the photon field in two dimensions has nontrivial topolog- 
ical properties formally expressed by the first homotopy 
group 7Ti[f7(l)] ~ Z. It implies the degeneracy of the 
ground state when each topological sector |n) is classified 
by integer number. Precisely this feature leads to non 
vanishing topological susceptibility with a "wrong" sign 
in two dimensions (5). The same degeneracy leads to 
nontrivial instanton solutions (8) interpolating between 
different topological sectors which saturate the topological 
susceptibility (13) with "wrong sign". 

In contrast to two dimensional case, in four dimensions 
one should not expect any contact term with a "wrong 
sign" similar to (5) as the third homotopy group is trivial, 
7T3(C/(1)) ~ 1, there is no degeneracy of the ground state 



as there is only a single trivial vacuum state. Therefore, 
one should not expect any non trivial surface terms in en- 
tropy computations in four dimensions. This expectation 
based on our conjecture is supported by explicit compu- 
tations [9] where it was shown that in four dimensional 
QED the surface term is gauge dependent and must be 
consistently discarded. In fact, we consider these argu- 
ments as further support for our conjecture formulated 
above as the behaviour of the system follows precisely the 
pattern dictated by this conjecture 8 . 

Our final comment is on interpretation of the surface 
term with a "wrong" sign given in conclusion of ref . [9] , 
where it has been suggested that, quote, "effective low- 
energy string theory which does not coincide with the 
ordinary QFT" in principle may produce some surface 
term with a "wrong" sign. We want to comment here, 
that in fact, very ordinary QFT may produce such kind of 
terms, which however, are non-dispersive in nature, and 
not related to any physical propagating degrees of freedom 
as explained in previous section II B for two dimensional 
case. As we argue in next section, such a behaviour is not 
a specific feature of two dimensional physics, but in fact 
very generic property in four dimensions as well. However, 
these nontrivial properties emerge in four dimensions only 
for non abelian gauge fields when the third homotopy 
group is non-trivial, ir 3 {SU(N)] ~ Z, the ground state is 
degenerate and each topological sector \n) is classified by 
integer number similar to two dimensional case considered 
in section II. The contact term with a "wrong sign" similar 
to eq. (5) is expected to emerge in this case as a result 
of nontrivial topological features of four dimensional non- 
abelian gauge theories. 



B. Four dimensional non-abelian QCD 

The goal of this section is to argue that all key elements 
from previous section II are also present in four dimen- 
sional QCD. In fact, the presence of the contact term 
with "wrong sign" in topological susceptibility in QCD 
is a crucial element of resolution of the so-called U(1)a 
problem [10, 11]. The difference with two dimensional 
case is that in strongly coupled QCD we can not perform 
exact analytical computations similar to (5,13). However, 
one can use an effective description in terms of auxiliary 
ghost field [11] to compute the non-dispersive contribution 
to topological susceptibility with "wrong sign" . This com- 



8 The argument presented above is based on observation that 4D 
space time (where computations [9] have been performed) has 
trivial topological properties. One can consider, instead, less 
trivial case when AD space time is represented, for example, by a 
torus, in which case the relevant homotopy group could be non- 
trivial, 7ri([/(l)) ~ Z, and contact term with a "wrong sign" in 
entropy computations may occur. In principle, this is a testable 
proposal. Technically, though, it could be quite a challenging 
problem. 



putation which employes the Veneziano ghost 9 is direct 
analog of derivation of eq. (19) when the Kogut Susskind 
ghost was used. Essentially, our goal here is to point out 
that the relevant features in 2d QED (discussed in section 
II when all computations can be explicitly performed) 
and in 4d QCD (when the final word is expected to come 
from the lattice numerical computations) are almost iden- 
tical. To further support these similarities we present 
some QCD lattice numerical results explicitly measuring 
the term with a "wrong sign" in topological susceptibility 
similar to eqs. (4,19). Based on these observations, our 
conjecture essentially implies that the entropy computa- 
tions in four dimensional non-abelian gauge theories must 
reveal a contribution with "wrong sign" as the crucial 
element, the degeneracy of the ground state, is present in 
the system. Moreover, it must be gauge invariant (and 
therefore, physical) in contrast with 4D QED computa- 
tions where it was shown to be gauge variant [9], and 
therefore, was discarded. 

Our starting remark is that the expression for topologi- 
cal density operator 

q = = ^z^ P *G a ^G apa = (28) 

being represented in terms of auxiliary scalar field $ has 
exactly the same form as in 2d Schwinger model, see 
section II C. The $ field in formula (28) is defined as 
K„ = and is the direct analog of representation (14) 
for 2d model. Our next remark is that four-derivative op- 
erator J d 4 x q 2 ~ / d 4 x (D$) 2 is expected to be induced 
in effective low energy Lagrangian as argued by Veneziano 
[11, 12] in his resolution of the U(1)a problem. As a 
result of generating of q 2 operator the relevant structure 
which emerges in the effective Lagrangian and describ- 
ing this system is identical to 2d QED case, i.e. it has 
precisely the same structure ~ $ [□□ + m^,D] <&. The 
corresponding path integral J 2?$ can be treated exactly 
in the same way as it was treated in 2d QED, i.e. it can 
be represented as the combination of the ghost <pi and a 
massive physical <j) field using the same trick by writing 
the inverse operator as follows 



aa + m 2 ,a m 2 ,\-n-m 2 , -u) ' y J 

in complete analogy with 2d case, see eq. (18). In fact, 
one can show that the relevant part of the low energy QCD 
Lagrangian in large N c limit in the form suggested by 
Veneziano [11, 12] is identical to that proposed by Kogut 
and Susskind for the 2d Schwinger model (15), where 
one should replace 2t/tt f^, 1 and fi —> m v i such that 

the scalar fields (j>i,<j>2,4> have appropriate (for 4D case) 
canonical dimension one, see [21] for details. This formal 
similarity leads to almost identical computations (in terms 



9 not to be confused with conventional Fadeev Popov ghosts. 
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FIG. 1. The density of the topological susceptibility x( r ) ~ 
(q(r),q(0)) as function of separation r such that x = j drx(r), 
adapted from [23] . Plot explicitly shows the presence of the 
contact term with the "wrong sign" (narrow peak around 
r ~ 0) represented by the Veneziano ghost in our framework. 



of the ghost) of the topological susceptibility in 4d QCD 
and in 2d QED. Indeed, by repeating all our previous steps 
leading to eq. (19) with known Green's functions which 
follow from (15) and with known expression for topological 
density operator q(x) ~ (\3<fi — 0<fi±) one arrives to the 
following expression for topological susceptibility 10 in 4d 
QCD in the chiral limit m = 0, 

XQCD = J d 4 x(T{q(x),q(0)}) (30) 

= ■ j d A x [5\x) - m 2 v ,D c ( mri ,x)] , 

where D c (m rj ix) is the Green's function of a 
free massive particle with standard normalization 
/ d 4 xm 2 , D c (m v > x) = 1. In this expression the S 4 (x) rep- 
resents the ghost contribution while the term proportional 
to D c (m v 'x) represents the physical rf contribution, see 
[21, 22] for details. The ghost's contribution can be also 
thought as the Witten's contact term [10] with "wrong 
sign" which is not related to any propagating degrees of 
freedom. The topological susceptibility XQCD( m = 0) = 
vanishes in the chiral limit as a result of exact cancellation 
of two terms entering (30) in complete accordance with 
WI. When m 7^ the cancellation is not complete and 
Xqcd — m (qq) a s it should. 



Of course x = to any order in perturbation theory because 
q(x) is a total divergence q = d^K^, However, as we learnt 
from [10, 11], x^0 due to the non-perturbative infrared physics. 
One can interpret field as a unique collective mode of the 
original gluon fields. It describes the dynamics of the degenerate 
states \n) representing the topologically nontrivial sectors of the 
ground state, it leads to a pole in unphysical subspace in the 
infrared, and finally, it saturates the contact term with a "wrong 
sign" in topological susceptibility (30). 



Similar to eq. (26) for two dimensional system, the 
non-dispersive term with "wrong sign" in topological 
susceptibility (30) in four dimensional QCD can be also 
represented as a surface integral 

X ~ / d 4 x[6 4 (x)] = J d^,(^). (31) 

In case of 2d QED we could compare our ghost's based 
computations (19) with exact results (4,5) and with ex- 
plicit summation over different topological sectors in pure 
EhM when no propagating degrees of freedom are present 
in the system (13). We do not have such a luxury in case 
of 4d QCD. Nevertheless, we can compare the ghost's 
based computations in 4d QCD given by eq.(30) with the 
lattice results, see e.g. [23]. We reproduce Fig.l from 
ref . [23] to illustrate few elements which are crucial for this 
work and which are explicitly present on the plot. First 
of all, there is a narrow peak around r ~ with a "wrong 
sign" . Second, one can observe a smooth behaviour in 
extended region of r ~ fm with the opposite sign. Both 
these elements are present in the lattice computations as 
one can see from Fig.l. The same important elements are 
also present in our ghost's based computations given by eq. 
(30). In different words, the QCD ghost does model the 
crucial property of the topological susceptibility related 
to summation over topological classes in gauge theories. 
This feature can not be accommodated by any physical 
asymptotic states as it is related to non-dispersive con- 
tribution in the topological susceptibility as explained 
above in section II B, and elaborated further in Appendix 
A where this feature is explained as a result of differences 
in definition between the Dyson's T-product and Wick's 
T- product. 

Our next step is to describe the behaviour of the 
same system (more precisely, the behaviour of the non- 
dispersive term in eq. (30, 31) proportional to the 5 i (x) 
function) in Rindler space in the presence of the horizon. 
We consider a simple case when the acceleration is suffi- 
ciently large Aq C£) 3> a 4 3> m\(qq) \ such that interaction 
term in (15) can be neglected and the Bogolubov's coef- 
ficients are exactly known. In this limit one can repeat 
all previous steps to arrive to the same Planck spectrum 
(25) for number density fluctuations of "fictitious parti- 
cles" with "wrong" commutation relations [18, 24]. This 
formula (up to some irrelevant numerical coefficient) has 
been reproduced in ref. [25] using a different technique. 
We interpret these fluctuations precisely in the same way 
as we did in section II D in 2d case when we interpreted 
these fluctuations of "fictitious particles" with "wrong" 
commutation relations as a different manifestation of the 
same physics which led to a "wrong sign" in entropy 
computations [5, 9]. As we emphasized before the cor- 
responding contribution is not related to any physical 
propagating degrees of freedom but rather, is due to topo- 
logical sectors in gauge theories which eventually lead to 
a non-dispersive contribution in topological susceptibility. 

Analysis of 2d case led us to the conjecture formulated 
at the end of section II that both phenomena ( "wrong sign" 
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in entropy computations and "wrong sign" in topological 
susceptibility) are originated from the same physics de- 
termined by surface dynamics of the "fictitious particles" . 
In this section we demonstrated that all relevant features 
are also present in 4d non-abelian QCD. Therefore, based 
on series of arguments presented above, it is naturally 
to assume that there will be a mismatch between black 
hole entropy and the entropy of entanglement in 4d QCD 
with the same "weird" features as listed in Introduction. 
However, in contrast with 4d QED we expect that the 
surface term with "wrong sign" in entropy will be a gauge 
invariant quantity similar to 2d QED case discussed in sec- 
tion II. This is essentially a prediction which follows from 
the conjecture. As we already mentioned there are other 
known examples when the degeneracy of the ground state 
in the presence of the horizon leads to mismatch between 
black hole entropy and entropy of entanglement, see Ap- 
pendix A for the details and references. One could argue 
that the dynamics of "fictitious particles" on a surface 
should be governed by the corresponding Chern-Simons 
action. We leave this subject for a future study [26]. 

IV. CONTACT INTERACTION AND 
PROFOUND CONSEQUENCES FOR 
EXPANDING UNIVERSE 

This portion of the paper is much more speculative 
in nature than the previous sections. However, these 
speculations may have some profound consequences on 
our understanding of expanding universe we live in when 
the horizon is inherent part of the system. Therefore, we 
opt to present these speculation in the present work. 

Non-dispersive contribution with a "wrong sign" in 
topological susceptibility (30) obviously implies, as eq. 
(1) states, that there is also some energy related to this 
contact term determined by the surface dynamics of "ficti- 
tious particles" . This 9— dependent portion of the energy, 
not related to any physical propagating degrees of free- 
dom, is well established phenomenon and tested on the 
lattice; it is not part of the debates. What is the part 
of the debates and speculations is the question on how 
this energy changes when background varies. In differ- 
ent words, the question we address in this section can 
be formulated as follows. How does the non-dispersive 
contribution to the 9— dependent portion of the energy 
vary when conventional Minkowski background is replaced 
by expanding universe with the horizon size L ~ H^ 1 
determined by the Hubble constant H 11 ? 



11 Here and in what follows we use parameter H ~ L~ 1 as a typical 
dimensional factor characterizing the visible size of our universe. 
We do not assume at this point that it is described by FRW metric 
with a single parameter H . In fact, it could be much more generic 
constructions when the spatial hypersurfaces are embedded in a 
compact 3d manifold such as, for example, Bianchi I geometry 
with few additional parameters. We refer to Appendix B of ref. 
[21] for a short review on this subject in a given context. 



The motivation for this question is as follows. We adopt 
the paradigm that the relevant definition of the energy 
which enters the Einstein equations is the difference AE = 
(E — EMink), similar to the well known Casimir effect when 
the observed energy is in fact a difference between the 
energy computed for a system with conducting boundaries 
(positioned at finite distance d) and infinite Minkowski 
space. In this framework it is quite natural to define the 
"renormalized vacuum energy" to be zero in Minkowski 
vacuum wherein the Einstein equations are automatically 
satisfied as the Ricci tensor identically vanishes. From 
this definition it is quite obvious that the "renormalized 
energy density" must be proportional to the deviation 
from Minkowski space-time geometry. This is in fact the 
standard subtraction procedure which is normally used 
for description the horizon's thermodynamics [27, 28] as 
well as in a course of computations of different Green's 
function in a curved background by subtracting infinities 
originated from the flat space [29] . In the present context 
such a definition AE = (E—EMmk) for the vacuum energy 
for the first time was advocated in 1967 by Zeldovich [30] 
who argued that p vac ~ Gm^ with m p being the proton's 
mass. Later on such a definition for the relevant energy 
AE = (E— -Eiviink) which enters the Einstein equations has 
been advocated from different perspectives in a number 
of papers, see e.g. relatively recent works [31-38] and 
references therein. 

This is exactly the motivation for question formulated 
in the previous paragraph: how does AE scale with HI 
The difference AE must obviously vanish when H — > 
as it corresponds to the transition to flat Minkowski 
space. How does it vanish? A naive expectation based 
on common sense suggests that AE ~ exp(— Aqcd/-^) ~ 
cxp(— 10 41 ) as QCD has a mass- gap ~ Aq CD , and there- 
fore, AE must not be sensitive to size of our universe 
L ~ H~ l . Such a naive expectation formally follows from 
the dispersion relations similar to (6) which dictate that 
a sensitivity to very large distances must be exponentially 
suppressed when the mass gap is present in the system. 

However, as we discussed at length in this paper, along 
with conventional dispersive contribution we also have 
the non- dispersive contribution (30, 31) which emerges as 
a result of topologically nontrivial sectors in four dimen- 
sional QCD. This contact term may lead to a power like 
scaling AE ~ H + 0(H 2 ) + ... rather than exponential 
like AE ~ exp(— Aqqd/H) because this term (in our 
framework) is described by massless ghost field (29) as 
discussed in section IIIB. The position of this unphys- 
ical massless pole is topologically protected as eq. (28) 
states, which eventually may result in power like scaling 
AE ~ H + 0(H 2 ) + ... rather than exponential like 12 . 



12 If a system is characterized by a single parameter, the curvature, 
then one should expect, on dimensional ground, that the first non- 
vanishing term in this expansion should be quadratic AE ~ H 2 
rather than linear. However, in a generic case one expects a linear 
non- vanishing term AE ~ H + 0(H 2 ) + see Appendix B of 
ref. [21] for the details. 
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In fact it was precisely the assumption postulated in [21, 
22] that the observable Dark Energy (DE) being identified 
with AE could be small but not exponentially small. 
Similar assumption based on very different arguments was 
also advocated in [31-34]. This postulate on Casimir- 
like scaling AE ~ H + 0(H) 2 has recently received a solid 
theoretical support as reviewed below. It is important 
to emphasize that this term with power like behaviour 
emerges as a result of non-dispersive nature of topological 
susceptibility (30), such that no violation of unitarity, 
gauge invariance or causality occur when the theory is 
formulated in terms of the unphysical ghosts [18]. If true, 
the difference between two metrics (expanding universe 
and Minkowski space-time) would lead to an estimate 

AE ~ HA 3 Q CD ~ (l(T 3 eV) 4 , (32) 

which is amazingly close to the observed DE value today. 
It is interesting to note that expression (32) reduces to 
Zeldovich's formula p vac ~ Gm® if one replaces Aqc d — > 
nip and H — > GAq CD . The last step follows from the 
solution of the Friedman equation 

H 2 = ^G ^ + ^ ^ pe>e ~ HA 3 Q CD (33) 

when the DE component dominates the matter compo- 
nent, Pde 3> Phi- In this case the evolution of the uni- 
verse approaches a de-Sitter state with constant expansion 
rate H ~ GAq CD as follows from (33). 

There are a number of arguments supporting the power 
like behaviour AE ~ H + O(H) 2 in gauge theories. First 
of all, it is an explicit computation in exactly solvable two- 
dimensional QED discussed in section II and defined in a 
box size L. The model has all elements crucial for present 
work: non-dispersive contact term (5) which emerges due 
to the topological sectors of the theory (13), and which 
can be described using auxiliary fictitious ghost fields 
(19). This model is known to be a theory of a single 
physical massive field. Still, one can explicitly compute 
AE ~ L^ 1 which is in drastic contrast with naively ex- 
pected exponential suppression, AE ~ e~ L [19]. It is 
important to emphasize that this correction AE ~ L^ 1 
while computed in terms of the ghost's (unphysical) de- 
grees of freedom in our framework, nevertheless represents 
a gauge invariant physical result. In different words, the 
final result AE ~ L^ 1 is not related to any violation of 
gauge invariance though it is computed using auxiliary 
fictitious ghost fields similar to computation of the contact 
term (19). 

One more support in power like behaviour is an explicit 
computation in a simple case of Rindler space-time in four 
dimensional QCD in the limit when a Rindler observer 
is moving with acceleration Aq CD 3> a 4 3> m\(qq)\ when 
the interaction term in eq. (15) can be neglected [18, 
24, 25]. These computations explicitly show that the 
power like behaviour emerges in four dimensional gauge 
systems in spite of the fact that the physical spectrum 
is gapped. In different words, a power like behaviour is 



not a specific feature of two dimensional physics as some 
people (wrongly) interpret the results of ref. [19]. 

Another argument supporting the power like corrections 
is the computation of the contact term in four dimensional 
QCD defined in a box size L. The computations are 
performed using the so-called instanton liquid model [39] . 
While the motivation for analysis [39] was quite different 
from our motivation, these model-based computations 
nevertheless explicitly show the emerges of power like 
corrections to non-dispersive portion of the topological 
susceptibility. 

Power like behaviour AE ~ L^ 1 is also supported by 
recent lattice results [40], see also earlier paper [41] with 
some hints on power like scaling in drastic contrast with 
naive expectations AE ~ exp — (Aqcd-^)- The approach 
advocated in ref. [40] is based on physical Coulomb gauge 
when nontrivial topological structure of the gauge fields is 
represented by the so-called Gribov copies. It is very dif- 
ferent from our approach where we advocate the auxiliary 
ghost's description to account for this physics. Eventu- 
ally, the physical results must not depend on the different 
technical tools which are used in different frameworks. 
However, it is not a simple task to demonstrate an inde- 
pendence of the results from an employed technique in a 
strongly coupled gauge theory! 

Finally, Casimir like scaling AE ~ L^ 1 can be tested in 
the so-called "deformed QCD" in weakly coupled regime 
when all computations are under complete theoretical 
control [42]. One can explicitly demonstrate that for the 
system defined on a manifold size L the 6- dependent 
portion of the energy shows the Casimir-like scaling E = 
—A ■ [l + j; + O (jji)] despite the presence of a mass 
gap in the system, in contrast with naive expectation 
E = —A'[l + B exp(— L)} which would normally originate 
from any physical massive propagating degrees of freedom 
consequent to conventional dispersion relations. 

Another remark worth to be mentioned is that the sign 
of AE = (E — i?Mink) is always expected to be negative 
in conventional quantum field theory computations. This 
is due to the fact that some modes can not be accommo- 
dated in a system with a nontrivial geometry/boundaries, 
and therefore the absolute value of -Etviink > E which cor- 
responds to AE < 0. The Casimir effect is the well known 
example when the sign (— ) emerges as a result of this 
subtraction procedure. The non-dispersive contribution 
into the energy, on the other hand, being represented by 
the ghost in our framework will lead to an opposite sign 
AE > 0. The positive sign for AE > is supported by an 
explicit computations in a simplified settings [18, 25, 42], 
and is consistent with observations corresponding to the 
accelerating universe with AE > 0. 

To conclude this section. The contact term with "wrong 
sign" in topological susceptibility which is present in gauge 
theories as a result of a nontrivial topological structure 
of the theory and which ultimately related to the "wrong 
sign" contribution in the entropy computations [5, 9], as 
argued in this paper, has another profound consequence. 
Namely, the very same physics, and the very same gauge 
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configurations may lead to a power like sensitivity from 
the distant regions AE ~ L^ 1 in drastic contrast with 
naive expectations AE ~ exp —(AqcD-k) which should 
occur from any conventional physical massive propagating 
degrees of freedom. If true, one can interpret the extra 
contribution to the energy (32), which we identify with 
DE, as a result of contact interaction with the horizon. 
This interpretation is consistent with interpretation of the 
term with a "wrong sign" in entropy computation in two 
dimensions [5, 9], and it is also consistent with our inter- 
pretation presented in section II D, see comments after 
formula (25). This interpretation is also consistent with 
arguments [5] suggesting that this term corresponds to a 
contact interaction with horizon in the description of the 
black hole entropy within a string theory formulation [43] . 



CONCLUSION. FUTURE DIRECTIONS. 

The main result of this paper is presented in form 
of a conjecture formulated at the end of section II and 
elaborated in section III. Essentially, the basic idea is 
that the surface term with a "wrong sign" in entropy 
computations [5, 9] and the contact term with "wrong 
sign" in topological susceptibility are both originated 
form the same physics, and both related to the same 
gauge configurations related to the nontrivial topological 
structure of the theory. If this conjecture turns out to 
be correct, it would unambiguously identify the nature of 
the well known mismatch between computations of the 
black hole entropy and entropy of entanglement for vector 
gauge fields. Similar mismatch (but in quite different 
context) was discussed also in [44-46]. In both cases the 
mismatch is a result of degeneracy of the ground state 
in the presence of the horizon, see Appendix A for the 
details. 

Another, much more profound consequence is that the 
same physics which is responsible for "wrong sign" contact 
term in topological susceptibility, and "wrong sign" con- 
tribution in the entropy computations [5, 9], may in fact 
lead to extra vacuum energy (32) (which is identified with 
observed DE) in expanding universe in comparison with 
Minkowski space. This extra energy emerges in gauge 
theories with multiple topological sectors as a result of 
merely existence of a causal horizon at distance L ~ H . 
Similar phenomenon as we already mentioned occurs in 
different systems [45, 46] when extra energy emerges as a 
result of dynamics of the "soft modes" at the horizon. The 
degeneracy of the vacuum state in the system discussed 
in [45, 46] is achieved by non- minimally coupling with a 
scalar field, see Appendix A for the details, while in our 
case the presence of topologically distinct sectors in the 
system is an inherent feature of the QCD dynamics. 

Here are some features of these unique gauge config- 
urations which are responsible for "wrong sign" in the 
entropy computations and "wrong sign" in topological 
susceptibility and which are characterized by very exotic 
properties drastically different from everything previously 



known: 

a) a typical wavelength of fluctuations of the auxiliary 
"fictitious particles" is determined by the horizon scale, 
Xk ~ 1/H ~ 10 Gyr, while smaller <C 1/H are ex- 
ponentially suppressed (25). Therefore, these modes do 
not gravitationally clump on distances smaller than the 
Hubble length, in contrast with all other types of matter, 
and can be identified with the observed properties of DE. 
Such very large wavelengths prevent us from adopting a 
meaningful scattering-based description, as the notion of 
particle is not even defined; 

b) the corresponding fluctuations are observer dependent, 
similar to the Unruh radiation, in contrast with any other 
types of radiation, see section II D and also [18] for detail 
discussions on problem of measurements in these circum- 
stances; 

c) the co-existence of the two drastically different scales 
(Aqcd ~ 100 MeV and H ~ 10~ 33 eV) is a direct conse- 
quence of the auxiliary conditions (17,23) on the physical 
Hilbert space rather than an ad hoc built-in feature such as 
small coupling or/and extra symmetries in a Lagrangian. 

So, essentially our proposal for the DE can be formu- 
lated as follows. The source for both: DE and mismatch 
between the black hole entropy and entropy of entangle- 
ment is the same and related to the dynamics of topo- 
logical sectors of a gauge theory in the presence of the 
horizon. In other words, the relevant gauge configura- 
tions which are responsible for DE are exactly the same 
which are responsible for the "wrong sign" -contribution 
in computations of refs [5, 9]. Precisely this contribution 
represents the mismatch between the black hole entropy 
and entropy of entanglement. In both cases the source 
for the extra term is the degeneracy of the vacuum state 
which is represented by different topological sectors, and 
in our framework is described by the Veneziano ghost. 
One should also add that a phenomenological analysis 
of the DE model based on this idea and represented by 
eqs.(32, 33) has been recently performed in [47] with 
conclusion that this model is consistent with all presently 
available observational data. 

Another important result of this work can be formulated 
as follows. When a system is promoted from Minkowski 
space to a an expanding universe, we expect a power like 
corrections AE ~ H + O(H) 2 rather than exponentially 
suppressed corrections, AE ~ exp — (Aqcd/^)- This 
happens in spite of the fact that the physical Hilbert 
subspace contains only massive propagating degree of 
freedom, and naively the sensitivity to very distant re- 
gions should be exponentially suppressed. However, the 
presence of the non-dispersive contributions (originated 
from degenerate topological sectors of the theory), which 
can not be associated with any physical asymptotic states 
falsifies this naive argument. Explicit computations in 2d 
QED [19], in 4d weakly- coupled "deformed QCD" [42] 
and numerical studies in real four dimensional QCD [40] 
where power like behaviour ~ L^ 1 is indeed observed, 
supports our claim. 

What is more remarkable is the fact that some of funda- 



13 



mental properties of gauge theories discussed in this paper 
can be, in principle, experimentally tested in relativistic 
heavy ion collider (RHIC) at Brookhaven and heavy ion 
program at LHC. In the "little bang" at RHIC the hori- 
zon appears as a result of induced acceleration a ~ Aqcd 
which itself emerges as a consequence of high energy col- 
lision. The acceleration a ~ Aqcd is a universal number 
which is determined by strong QCD dynamics, does not 
depend on energy or other properties of the colliding par- 
ticles, and plays the role of Hubble constant H ~ 10~ 33 eV 
of expanding universe, see [24] for the details. 
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Appendix A: Subtleties in definition of energy. 

It is quite natural to expect that there should be some 
energy associated with fluctuations of fictitious particles 
(25). As we argued above these fluctuations, which reflect 
the nontrivial topological structure of the gauge theory, is 
the source for missmatch between the black hole entropy 
and entropy of entanglement in the presence of the horizon. 
It is clear that there is some ambiguity in definition of 
this type of energy due to a number of reasons. First 
of all, as we discussed in the text the relevant physics 
is determined by the dynamics on the surface rather 
than in the bulk of the space-time. Therefore, it is not 
obvious if a simple insertion Uk into definition (22) would 
properly reflect this feature. Another comment is that 
the physics of fluctuations (25) is the observer dependent 
property similar to the Unruh effect as discussed in section 
II D. Therefore, all subtleties related to the Unruh effect 
are also present here. Also, very large wavelengths of 
these fluctuations prevent us from adopting a conventional 
description in terms of particles, as the notion of particle 
is not even defined. 

Finally, and most importantly, the contact interaction 
which is the main subject of this paper, can not be ex- 
pressed in terms of any physical states, but rather is 
formulated in terms of fluctuations of fictitious particles 
with "wrong" commutation relations. These unphysical 
states contribute to the non-dispersive portion of the 
correlation function, not to the dispersive part which is 
unambiguously determined by physical spectral function 
through conventional dispersion relations. Our descrip- 
tion in terms of the ghost is simply a convenient way to 
study this IR physics in a covariant gauge. The same 
physics in Coulomb (physical) gauge when the ghost de- 



grees of freedom are not present in the system leads to the 
long range forces as discussed in simple two dimensional 
model long ago [16]. These long range forces prevent 
distant regions from acting independently. The vacuum 
in this system is degenerate, and one should consider 
an infinite superposition of of the winding states \n) as 
originally discussed in [16]. We think that precisely this 
feature prevented the author of ref. [5] to use the phys- 
ical Coulomb gauge in two dimensions in computations 
of the entropy, where a covariant gauge has been used 
instead. The same physics in Coulomb gauge in 4d QCD 
is reflected by existence of the so-called Gribov copies, 
and one should use some numerical lattice methods to 
study the relevant physics [40]. 

The reason why we pay so much attention to the topo- 
logical susceptibility x and the corresponding contact 
term ~ S(x) which enters the expressions for topological 
susceptibility (4) and (30) is due to its relation to the 9 
dependent portion of the vacuum energy as eq. (1) states. 
Therefore, the presence of non-dispersive contributions 
in x automatically implies presence of the corresponding 
non-dispersive contribution in the vacuum energy E vac . 
At the same time as we discussed in the text the non- 
dispersive contribution in x (and in the vacuum energy) 
can be interpreted as a result of nontrivial topological 
structure of the gauge theories. These discussions explic- 
itly demonstrate some subtleties on possible definition of 
the vacuum energy which should accommodate the physics 
related to the contact term discussed in this paper. 

It is quite fortunate that in specific case with compu- 
tations of x we can easily separate non-dispersive and 
physical contributions. This is due to the fact, that these 
two very different terms contribute to x with opposite 
signs as discussed in the text. Indeed, the topological 
susceptibility in pure YM gauge theory (no quarks, and 
no rf contribution) according to [10, 11] is given by 13 

Xym =i lim / d*xe ikx (T{q(x),q(0)}) w = 

1 , Tin 1 

= -X 2 YM <0, A^ M = J Y L , (Al) 

where (...)w stands for Wick T- product, see below. The 
expression on r.h.s of eq. (Al) corresponds to the sub- 
traction constant ~ S 4 (x) in eq. (30). This term has 
a "wrong sign" (in comparison with contribution from 
any real physical states), the property which motivated 
the term "Veneziano ghost" . Indeed, a physical state of 
mass ma, momentum k — >■ and coupling (0|g|G) = cq 
contributes to the topological susceptibility with the sign 



13 All formulae in this Appendix A are written in Minkowski space in 
contrast with our discussions in the text, where a comparison with 
path integral computations and lattice numerical computations 
(which arc always performed in Euclidean space) was made. In 
particular, there is factor "i" in the definition of the correlation 
function (Al). 
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which is opposite to (Al), 

ilim [ d 4 xe ikx (T{q(x),q(0)}) D ~ (A2) 

fc-s-0 j 

ilim(0|g|G) , 2 1 2 (G\q\0) ~ J^j! > 0, 

where (...)d stands for Dyson T- product, see below. How- 
ever, the negative sign for the topological susceptibility 
(Al) is what is required to extract the physical mass for 
the 77' meson, see the original reference [12] for a thor- 
ough discussion. The difference between behaviour (Al) 
and (A2) is related to inequivalent definitions of these 
correlation functions. The behaviour (A2) corresponds to 
the usual Dyson T-product when only physical states con- 
tribute, while eq. (Al) corresponds to Wick T-product 
obtained by variation of the partition function over 9 
parameter. The difference in the definitions constitutes 
precisely the subtraction constant ~ S i (x) in eq. (30). 
The WI expressed as XQCD(m q = 0) = is satisfied for 
the Wick T-product, not for the Dyson T-product. 

It is interesting to note that an analogous phenomenon 
(but in quite different context) was discussed in ref. [45, 
46], where it was observed that there are two definitions 
of energy when the difference is saturated by the so-called 
"soft modes" fluctuating in far infrared at the horizon. 
First definition is the canonical energy determined by 
the hamiltonian which is generator of translations of the 
system along the timelike Killing vector field £ M . Since 
the Killing vector £ M vanishes at the bifurcation surface 
of the Killing horizons, the corresponding Hamiltonian 
is degenerate. Therefore, one can add to the system an 
arbitrary number of "soft modes" without changing the 
canonical energy. These "soft modes" contribute to the 
surface integral which was interpreted as a Noether charge 
of some non-minimally coupled scalar field <p s . Precisely 
this contribution of the "soft modes" distinguishes two 



different definitions of the energy. As argued in ref. [45, 
46] precisely the dynamics of "soft modes" constitutes the 
difference between Bekenstein -Hawking entropy and the 
entropy of entanglement. 

It is very similar to our case when we argued that a 
"wrong sign" in entropy computations of refs [5, 9] is a 
result of degeneracy in the gauge theory represented in our 
framework by the fluctuations of the fictitious particles 
with a typical wavelength of order of the horizon scale, 
Xk ~ H~ l . In different words, the "soft modes" from 
[45, 46] play the same role as pure gauge fields which 
describe different topological sectors of the theory in our 
case. In both cases the effect emerges as a result of the 
degeneracy of the ground state in the presence of the 
horizon, and in both cases the difference is determined 
by some surface integrals. Further to this analogy, two 
different types of energies which can be reconstructed 
from two different definitions of topological susceptibility 
(Al) and (A2) precisely correspond to two different types 
of energies discussed in refs. [45, 46]. 

The final comment on the similarities between the two 
very different systems is as follows. The relevant gauge 
configurations which are responsible for mismatch be- 
tween the black hole entropy and entropy of entanglement 
from computations [5, 9] are exactly the same which are 
responsible for the contact term in topological suscepti- 
bility (30). This mismatch is very similar to extra term 
discussed in [45, 46] which is resulted from the dynam- 
ics of the "soft modes" at the horizon. In both cases 
the mismatch is a result of degeneracy present in both 
systems. In refs. [45, 46] this degeneracy is a result of 
non-minimally coupling with scalar field <p s . In our cases 
this degeneracy is a reflection of topological structure of 
gauge theories, and inherent feature of QCD. However, 
the outcome of the degeneracy is very similar in both 
cases as we argued above, and can be interpreted as a 
contact interaction with horizon, see discussions at the 
end of section IV. 
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